Microscopic origin of magnetism and magnetic interactions in ferropnictides 
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One year after their initial discovery, two schools of thought have crystallized regarding the 
electronic structure and magnetic properties of ferropnictide systems. One postulates that these are 
itinerant weakly correlated metallic systems that become magnetic by virtue of spin-Peierls type 
transition due to near-nesting between the hole and the electron Fermi surface pockets. The other 
argues these materials are strongly or at least moderately correlated, the electrons are considerably 
localized and close to a Mott-Hubbard transition, with the local magnetic moments interacting 
via short-range superexchange. In this paper we argue that neither picture is fully correct. The 
systems are moderately correlated, but with correlations driven by Hund's rule coupling rather than 
by the on-site Hubbard repulsion. The iron moments are largely local, driven by Hund's intra- 
atomic exchange. Superexchange is not operative and the interactions between the Fe moments are 
considerably long-range and driven mostly by one-electron energies of all occupied states. 

PACS numbers: Pacs Numbers: 



Ferropnictides are still attracting widespread attention 
from researchers both inside and outside the field of su- 
perconductivity. There is now a nearly universal agree- 
ment that magnetism and, specifically, proximity to an 
antiferromagnetic "stripe-order" transition plays a ma- 
jor role in the physics of these compounds. There is also 
growing evidence that the magnetic properties and corre- 
lation effects in this system are not controlled by the Hub- 
bard U as in cuprates (spectroscopy tells us that Hubbard 
correlations are weak and the effective U is on the order 
of 1 eVd, smaller than the bandwidth; first principles 
calculations of U support thisQ])- On the other hand, 
the multiband character of Fe bands and the large intra- 
atomic (Hund's) exchange coupling in Fe suggest that the 
Hund's J may play the main role in the magnetism.. 

As opposed to the Hubbard U, the Hund's J is gener- 
ally well accounted for in density-functional calculations 
(where it is called the Stoner /).) Indeed, the local den- 
sity approximation (LDA) correctly predicts the particu- 
lar antiferromagnetic and structural ground state of un- 
doped ferropnictides, in striking contrast to the cuprates. 
In view of this, it is instrumental to trace down the origin 
of magnetism within LDA, and to disentangle the na- 
ture of magnetic interactions captured by this approach. 
It is highly likely that the physics uncovered by density 
functional theory will reflect the actual physics of these 
systems. Given the heated (but largely devoid of solid 
facts) discussion of whether antiferromagnetism in pnic- 
tides is due to Fermi surface nesting or to second neighbor 
superexchange (See Ref. Q for a review), a clear under- 
standing of, at least, the message that LDA calculations 
send seems highly necessary. 

In this paper, we analyze the magnetic interactions and 
demonstrate that neither of the above two views (often 
presented as an axiomatic dilemma) is correct. The mag- 
netism appears due to local Hund's rule coupling, while 
the particular ground state is selected by itinerant, essen- 
tially one-electron interactions, of which only a small part 



is played by the Fermi surface nesting. We will also ex- 
plain why conventional "Anderson-Kanamori" superex- 
change is not operative here, and will show some striking 
examples where calculations and experiment contradict 
both local superexchange and spin-Peierls pictures, yet 
are perfectly understandable on the basis of one-electron 
energy balance. 

We start with a qualitative analysis. The Hund's rule 
coupling energy in density functional theory is expressed 
as Eh — —Im^/A, where m is the magnetization of 
an Fe ion, and / « 0.8 eV is the Stoner factor for Fe 
(0.9 eV in GGA). Depending on the material, the self- 
consistent magnetic moment on Fe appears to be between 
1.5-2 hb in LDA and 1.8-2.5 in GGA. The corresponding 
energy gain even in LDA is 0.5 eV, which is remark- 
ably large. In other words, every individual Fe wants 
to be strongly magnetic and the advantage of spin po- 
larization should lead to a magnetic ground state at the 
mean field level, unless an unusually large kinetic energy 
penalty exists. However, this is exactly the case for the 
formation of a ferromagnetic configuration. To create 
a magnetization m on Fe, one needs to move approx- 
imately 1.15 (to account for the relative share of Fe-d 
orbitals at the Fermi level) spin-minority electrons into 
unoccupied spin-majority states, incurring an energy loss 
of « (1.15m)ViV|(£;F). The density of states (DOS) per 
Fe, N^{Ep), varies between 1 and 1.5 eV~^, depending 
on the system, creating an energy loss for m = 1.5 /is 
of 0.5-0.8 eV. This cost is about as large as the Hund's 
rule energy gain estimated above. This shows that the 
system is on the verge of a ferromagnetic instability, but 
nothing more. 

In low-DOS metals, magnetization without a large cost 
in kinetic energy is possible if some type of antiferro- 
magnetic arrangement is formed (cf. metal Cr and Mn). 
For a broad band metal, this narrows the conductivity 
band, but as long as the exchange splitting is smaller 
than the bandwidth, the cost is small. Because in fer- 
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FIG. 1: (color online) Top-down view of the a) checkerboard 
b) stripe and c) doublestripe magnetic patterns for a single 
FeAs or FeTe layer. The light colored sites have majority up 
spin and the darker sites have majority down spin. 

ropnictides the calculated bandwidth is 5-6 eV and the 
exchange splitting ml is at most 2 eV, this mechanism 
should be very favorable. 

It is interesting to consider how the system determines 
which particular AFM arrangement is most profitable 
from the point of view of the one-electron energy (note 
that LDA calculations can be forced to converge to nearly 
any AFM pattern, but not to a ferromagnetic state). If 
the resulting magnetization is small, the answer is obvi- 
ous: the second derivative of the total energy with respect 
to magnetization is defined by the noninteracting suscep- 
tibility at the AF wave vector Q, d'^E/dvn? — — Xo ^(Q) 
(with the small caveat that an actual spin density wave 
is not a single harmonic, but includes all wave vectors 
Q+G, where G is a reciprocal lattice vector). The imag- 
inary part of xo is directly related to Fermi surface nest- 
ing, being defined, in the constant matrix elements ap- 
proximation, as X^ij / '5(^ki)'5(£k+Q,i)rfk, while the (ac- 
tually relevant) real part collects information from all 
states and may or may not have any relation to the nest- 
ing conditions (for a detailed discussion see Ref. Q). 

Geometrical nesting, as a property of the Fermi sur- 
face, becomes even more disconnected from a real insta- 
bility in the strongly nonlinear regime, ?n 1 /is, which 
is the case for ferropnictides. Monitoring the evolution of 
the electronic bands with increasing spin polarization[5], 
one observes that at m ~ 1 /is the resulting bands can 
in no way be described as anticrossing downfolded non- 
magnetic bands with partial gapping of the Fermi sur- 
face. Rather, the entire Fe d band is fully restructured. 
Although the lowest-energy AFM state wave vector in- 
deed coincides with the quasi-nesting wave vector in some 
cases, it is not always true, as exemplified by the case of 
FeTe that we discuss later. 

It should be noted that while quasi-ncsting is not par- 
ticularly relevant for the long-range ordering in the un- 
doped crystals, it does define the low-energy excitations 
in non-magnetic phases and these can perfectly well me- 
diate superconductivity. 

Having established a general framework, we now ad- 
dress specific examples. First, we investigate checker- 
board, stripe, and double-stripe magnetic structures (See 
Fig. [T]) and show that the stripe order is lower in en- 
ergy than either the checkerboard or the double-stripe 



TABLE I: Stabilization energies for various magnetic config- 
urations in the 122 and FeTe systems. All energies are per Fe 
atom. 





checkerboard 


stripe 


double stripe 


BaFe2As2 
FeTe 


16 meV 


94 meV 
207 meV 


0.6 meV 
230 meV 



structure for the 122 systems, but not for FeTe. We use 
BaFe2As2 as an example, but the results for LaFeAsO 
are very similar. Our calculations were performed using 
an all-electron, full-potential LAPW package WIEN2k, 
in the Generalized Gradient Approximation, similar to 
Ref. Q. All structures were fully relaxed (except where 
stated otherwise) using the Vienna Ab-initio Simulation 
Program (VASP) "9], with the PAW formulation ^10] and 
also using GGA. In Table 1 we show the magnetic stabi- 
lization energies of the three different antiferromagnetic 
structures. 

In Fig. [2^,b,c, we show the DOSs for BaFe2As2 in each 
of the three magnetic configurations along with the non- 
magnetic DOS. Compared to the nonmagnetic DOS, we 
see that the checkerboard pattern has a very similar spec- 
trum at and near the Fermi energy and gains one-electron 
energy by shifting spectral weight from the region be- 
tween -0.5 and -1.0 downward to the region between -1.0 
and -2.0. The doublestripe pattern actually incurs an 
energy penalty at and just below Ei?, but gains energy 
by shifting weight downward from between -0.2 and -0.7 
to between -1.0 and -2.0. The ground state configura- 
tion, in contrast to the other two, gains energy all the 
way from E^? to -0.9 by shifting weight downward. This 
is accomplished through the opening of a large pseudo- 
gap (this terminology has no connection with the pseu- 
dogap in cuprates and simply signifies a depression in 
one-electron DOS around the Fermi level). Though all 
three magnetic configurations are stable with respect to 
a nonmagnetic state, it is visibly the case that the stripe 
ordering has the greatest one electron energy advantage. 
This is reflected in the much larger gain in total energy 
(See Table U). 

Let us now compare the results with the same calcu- 
lations for FeTe. As indicated in a number of papers, 
FeTe is always formed with an excess Fe, so the fact 
that experiment gives the double stripe structure as the 
low-temperature ground state [11] should be taken cum 
grano salts. However, as Table 1 shows, it is definitely 
the stoichiometric ground state in density functional cal- 
culations, and this is the only thing that matters for our 
analysis [l^j ■ We note here that we do not fully relax the 
FeTe structure, but only relax the internal positions. As 
before these relaxations are done separately for magnetic 
and nonmagnetic cases. 

Table |T] indicates that for FeTe, as opposed to 
BaFe2As2, the energies of the single and double stripe 



3 







a) 


non-magnetic 

checkerboard 






b) 


non-magnetic 

stripe 






c) 


non-magnetic 

double stripe 


1 


1 . 



Energy (eV) 



FIG. 2: (color online) The densities of states for BaFe2As2 in 
the non-magnetic configuration in comparison to a) checker- 
board magnetic pattern b) stripe (ground state) magnetic pat- 
tern and c) double stripe magnetic pattern 
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FIG. 3: (color online) The densities of states for FeTe in the 
non-magnetic configuration in comparison to a) stripe mag- 
netic pattern and b) double stripe magnetic pattern. 

phases are relatively close. This suggests that the crude 
method of determining the ground state by looking at 
the DOS may not work here, as the DOSs for the two 
AFM phases will probably be similar. Indeed, this is 
what we see in Fig. [3K,b where a large downward shift of 
spectral weight is visible for both patterns. Interestingly, 
the nonmagnetic Fermi surfaces in FeTe are extremely 
similar to those in the 122 and 1111 materials, whereas 
nonmagnetic DOS and the magnetic electronic structure 
are quite different. This reinforces that FS nesting, which 
would be nearly identical for BaFe2As2 and FeTe, is not 
driving the magnetic order. 

Our calculations also provide a strong argument 
against superexchange. Looking at the patterns in Fig. 
[U it is easy to see that for both the stripe and double 



stripe patterns, the first neighbor exchange, Ji, does not 
contribute, due to equal numbers of aligned and anti- 
aligned spins. The second neighbor exchange, J2, would 
have to be stronger than Ji/2 in order for the stripe pat- 
tern to be energetically favorable over the checkerboard 
pattern in a superexchange picture. It has often been 
argued that this situation is not unreasonable since the 
Fe-As-Fe paths available for Ji and J2 are similar. For 
double stripe order in the FeTe system, however, both Ji 
and J2 cancel, leaving only J3 to establish the ordering. 
Considering the remarkably strong stability (compared 
to non-magnetic) calculated for double stripe order (See 
Table HJ, this is hard to rationalize. Furthermore, the 
energy term for stripe is J2 - J3 (compared to J3 alone 
for double stripe). For double stripe to stabihze, J3 could 
be no smaller than J2/2, but the "similar hopping paths" 
justification used for Ji and J2 and is not available: the 
third neighbor exchange path is more than twice as long 
as the second neighbor one and involves As- As hopping 
across a distance of a full lattice constant. Thus, the ex- 
istence and stability of the double stripe order severely 
strains the credibility of the superexchange picture. This 
is, in fact, to be expected since superexchange is not efh- 
cient when the bandwidth is much larger than the energy 
cost of flipping an electron's spin, which is precisely the 
case here. 

This does not, however, mean that one cannot map 
the dependence of the total energy onto a suitable short- 
range exchange model. In fact, it is hard to imagine a 
case in which this would not be possible. Yet, in carrying 
out this procedure for ferropnictide systems, one should 
be aware of the following caveats: 

(1) There is no microscopic justification (as for instance 
in the Hubbard model) for introducing any J — t (or 
Ji — J2 — t) Hamiltonian. 

(2) There is no guarantee that this kind of mapping 
can be stopped at first or second neighbors. In fact, accu- 
rate calculations show that at least some of the exchange 
parameters in these mappings decay as l/i?"^, just as in 
metal ironQ. 

(3) The resulting exchange parameters strongly and 
qualitatively depend on the long-range order established 
in the system. In particular, the parameters that can 
be used to describe the ordered state cannot be used to 
describe the spin fluctuations, and vice versa. (See Ref. 
1^ and references therein.) 

(4) In the absence of superexchange, there is no reason 
to believe that the total energy can be mapped onto a 
Heisenberg model, ^j^- JijSi ■ Sj. In fact, direct calcula- 
tions show that at least one biquadratic term needs to 
be added to map the total energy onto the mean-field 
Hamiltonian, ^2^(8^ ■ S.+i)^, where X - J. ^ 

We now switch our attention to the structural transi- 
tion observed simultaneoulsy with the magnetic one in 
the 122 systems. Density functional calculations very ac- 
curately reproduce the experimentally observed distor- 
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with respect to the magnetic ordering and structural dis- 
tortion in the ferropnictides hes in the one-electron en- 
ergies. Our results resolve the superficially binary choice 
between superexchange and Fermi surface nesting in fa- 
vor of a third mechanism that is neither fully local- 
ized nor fully itinerant. One-electron energy is gained 
throughout an energy range of at least 1 eV below E^- 
and the ground state is determined by which magnetic 
pattern most effectively exploits a downshift in spectral 
weight, not by fermiology. On the other hand, the Fermi 
surface itself is the operative feature for determination 
of the structural distortion. The energy minimum for an 
in-plane distortion corresponds to a simultaneous mini- 
mization of the Fermi surface area. 

We thank H. Eschrig, K. Koepernik and T. Yildirim for 
useful and engaging discussions related to this work. We 
acknowledge funding from the Office of Naval Research. 



FIG. 4: The Fermi surfaces of stripe-ordered BaFe2As2. Top 
panel shows the 'reverse distortion' in which the Fe-Fe dis- 
tance is lengthened along the like spin direction and shortened 
along the unlike spin direction. The bottom panels shows the 
fully relaxed calculation which reproduces the experimentally 
observed distortion (to within a few percent). 



tion in which Fe ions along the stripe direction are closer 
to one another than Fe atoms belong to adjacent stripes 
I3I . 14 1 . We investigated whether the structural distor- 
tion, like the magnetic ordering, can be understood in 
terms of one electron energies by calculating the DOS 
for a variety of small changes in the a and b lattice con- 
stants. In contrast to changing the magnetic pattern, 
changing the structural distortion has very little effect 
on the DOS away from the Fermi energy. There were no 
large shifts of spectral weight to lower energies, though 
small shifts of the order of 0.05 eV did occur and these 
were within 0.5 eV of the Fermi energy (for comparison 
see the heavy restructuring of the DOS in Figs. [2]and[3]). 
The distortion can therefore be treated as a linear pertur- 
bation with a one-electron energy lowering observable at 
(or very near) the Fermi energy. Specifically, we find that 
the lowest energy structure corresponds to the smallest 
Fermi surface area. As an example, in Fig. [4] we show 
the Fermi surface in the magnetic Brillouin zone of the 
fully relaxed (lowest energy) structure and a 'reverse dis- 
tortion' in which the distances between like and unlike 
spins are reversed from the correct configuration. The 
change in the size of the Fermi surface is clearly visible. 
We were unable to engineer a further minimization of 
the Fermi surface with any choice of in-plane distortions 
other than the optimal energy one. 

In conclusion, we have shown that the relevant physics 
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